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Different theoretical approaches to the famous two state Landau - Zener problem are briefly dis- 
cussed. Apart from traditional methods of the adiabatic perturbation theory, Born - Oppenheimer 
approximation with geometric phase effects, two-level approach, and momentum space representa- 
tion, the problem is treated semiclassically also in the coordinate space. Within the framework of 
the instanton approach we present a full and unified description of ID Landau-Zener problem of 
level crossing. The method enables us to treat accurately all four transition points (appearing at two 
levels crossing) , while the standard WKB approach takes into account only two of them. The latter 
approximation is adequate for calculating of the transition probability or for studying of scattering 
processes, however it does not work for finding corresponding chemical reactions rates, where very 
often for typical range of parameters all four transition points can be relevant. Applications of the 
method and of the results may concern the various systems in physics, chemistry and biology. 
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I. INTRODUCTION 

At first sight, the title of this paper might sound perplexing. What else can be said about Landau - Zener (LZ) 
problem after numerous descriptions in both research and textbook literature? However, although theoretical (and 
experimental) investigations of different LZ systems began more than seventy years ago, it is still remains an active 
area of research. Various approaches to LZ problem have appeared separately in the literature (see e.g. by no means 
not a complete list of publications [1] - [67]), these references are not fully consistent with each other. We, therefore, 
think that it is important to collect and discuss all these approaches in one place. We will study the ID LZ problem 
[1] of quantum mechanical transitions between the levels of a two-level system at the avoided level crossing. LZ 
theory treats a quantum system placed in a a slowly varying external field. Naturally in such a condition the system 
adiabatically follows variation of an initially prepared discrete state until its time dependent energy level crosses 
another one. Near the crossing point evidently the adiabaticity condition is violated (like a semiclassical behavior 
is violated near turning points). The slow variation of the perturbation means that the duration of the transition 
process is very long, and therefore the change in the action during this time is large. In this sense the LZ problem is 
a semiclassical one (but with respect to time instead of a coordinate for standard semiclassical problems). 

As is well known the problem presents the most basic model of non-adiabatic transition which plays very important 
role in many fields of physics, chemistry, and biology. Not surprisingly therefore that numerous monographs and 
uncountable papers have been devoted to this subject. In the literature there are roughly speaking three kinds of 
semiclassical modelling of LZ problem, namely: 

• (i) two level system approach [2] - [8]; 

• (ii) adiabatic perturbation theory [9] - [21] (see also review paper [6]) ; 

• (iii) momentum space representation [22] - [25] . 

Because different approaches have been proposed to study the LZ problem one of the immediate motivation of the 
present paper is to develop an uniform and systematic procedure for handling the problem. We will show that all 
three methods are equivalent for treating tunneling and over-barrier regions of parameters, and no one of them can 
be applied to study say intermediate region of parameters, where all four of the involving into LZ system states are 
relevant. To study this region is our main concern in this paper. A second question addressed here concerns so-called 
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connection matrices. At usual textbook treatment of the LZ problem, only the transition probabilities are calculated 
and expressed in terms of the genuine two-level LZ formula successively applied at each diabatic level intersection. 
Evidently such a procedure is an approximation for the general LZ problem including even in the simplest case at 
least 4 energy levels. To solve many important physical or chemical problems one must find the 4x4 (not only 2x2 
connection matrices related these 4 states. 

While our paper is not intended as a comprehensive review we detail here key results of the standard WKB and 
instanton approaches from our own researches and literatures within the context of different factors that we feel 
are important to studying LZ problem. Specifically we focus in the next section (II) on the Born - Oppenheimer 
approximation which is a benchmark to test semiclassical approximations. In section (III) we lay the foundation for 
treating LZ problem, namely - adiabatic perturbation theory. Section (IV) is devoted to the generalization of the 
instanton method enables us to investigate LZ problem in the momentum space. It is shown in this section that for a 
linear (in a ID coordinate under consideration) potential WKB semiclassical wave functions in the momentum space 
coincide with the instanton wave functions. For the quadratically approximated (parabolic) potentials the instanton 
wave functions arc exact and have no singularities (unlike WKB wave functions; remind that the same kind of relations 
hold for the WKB and instanton wave functions in the coordinate space [26] - [29]). 

In this paper we are advocating for the instanton approach, but it is worth noting that, nevertheless, many important 
results have been obtained in the frame work of the WKB approach [1] -[8]. For example, one of the very efficient 
technique (so-called propagator method) was proposed and elaborated by Miller and his coworkers [34], [35], [36] (see 
also [26]). This approach uses semiclassic (van Fleck - Gutzwiller types) propagators, taking into account automatically 
in terms of the general WKB formalism, the contribution coming from the contour around a complex turning point. 
The accuracy of the WKB method can be improved considerably, [30] [2], [31], [5] (more recent references on so-called 
Laplace contour integration can be found also in [32]) by the appropriate choice of the integration path around the 
turning point, and it appears to be quite accurate for the tunneling and over-barrier regions, but it becomes non- 
adequate in the intermediate energy region. It has been overlooked in the previous investigations treating this region 
by a simple interpolation from the tunneling (with monotonic decay of the transition probability) to the over-barrier 
(with oscillating behavior) regions. 

We present all details of the LZ problem for two electronic states in section V using the instanton description of 
LZ problem in the coordinate space. In the section the basic two second order differential (Schrodinger) equations 
to be dealt with are written in the so-called diabatic state representation (i.e. in the basis of "crossed" levels). 
Neglecting higher order space derivatives we find asymptotic solutions, and using adiabatic - diabatic transformation 
we match the solutions in the intermediate region. The complete scattering matrix for the LZ problem is derived 
in the section VI. In the section VII we derive the quantization rules for crossing diabatic potentials and discuss 
shortly the application of the obtained results to some particular models of level crossings which are relevant for 
the interpretation and description of experimental data on spectroscopy of non-rigid molecules, on inelastic atomic 
collisions [33], non-radiative transitions arising from " intersystem" crossings of potential energy surfaces in molecular 
spectroscopy and chemical dynamics (see e.g. [26] and references herein). In the last section VIII we draw our 
conclusions. 

Throughout what follows we will consider ID case only. The LZ problem for ID potentials coupled with the thermal 
bath of harmonic oscillators is shown to reduced to a certain renormalization of the Massey parameter, where entering 
the expression for this parameter longitudinal velocity is decreased due to coupling to transverse oscillations (see [26] 
and references herein, and for more recent references also [66], [67]). Of course the energetic profile of any real system 
is characterized by a multidimensional surface. However, it is often possible to idc;ntify a reaction coordinate, such that 
the energy barrier between initial and final states is minimized along this specific direction, and, therefore, effectively 
one can treat the system under consideration as ID. In certain systems, the physical interpretation of the reaction 
coordinate is immediate (e.g. the relative bond length in two diatomic molecules), but sometimes it is not an easy 
(if possible at all) task, because of the many possibilities involved. The latter (multidimensional) case will be studied 
elsewhere. Unfortunately the accuracy of the WKB method near the barrier top is very poor to make any numbers 
realistic and it is one more motivation to use the alternative to WKB semiclassical formalism - extreme tunneling 
trajectory or instanton technique. 

II. BORN - OPPENHEIMER APPROXIMATION. 

It may be useful to illustrate the essential physics of the LZ problem starting with a very well known picture 
corresponding to the Born - Oppenheimer approximation [37] , [1] which leads to the separation of nuclear and electronic 
motions, and the approximation is valid only because electrons are so much lighter than nuclei and therefore move so 
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much faster. Thus the small parameter of the Born - Oppenheimer approximation is 



A: 



/me\ 1/4 
— «1 

V m / 



(2.1) 



where me and m are electronic and nuclear masses respectively. On the other hands the semiclassical parameter 



7 



> 1 



(2.2) 



where a is a characteristic length of the problem, and the characteristic nuclear vibration frequency oc m~^/^, 
therefore 7 oc A~^. Prom this simple fact important conclusions are arrived at. Indeed one can satisfy the semiclassical 
condition 7 ^ 1 by assuming formally ?i ^ or equivalcntly A ^ 0. This correspondence allows us to apply on the 
same footing to the separation of scales for nuclear and electronic motions either the Born - Oppenheimer or the 
semiclassical approximation. 

In the traditional Born - Oppenheimer approach the solution ^' to the full (i.e. including electronic Hamiltonian 
ffg) depending on electronic coordinates r, and nuclear Hamiltonian depending on nuclear coordinates R) Schrodinger 
equation is presented in the form of an expansion over the electronic Hamiltonian eigen functions (f)n 



* = ^$„(i?)0„(r, R). 



(2.3) 



The electronic cigcn values En depends on the nuclear coordinates, and the expansion coefficients ^n{R) is determined 
by the Born - Oppenheimer equations 
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2m ^ 



where for to 7^ A; 
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^ {5nk^R-iAnk){hmyR-iAkm)4'm, (2.4) 
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•,k = n 



»|V_R(/>ft) , 



(2.5) 



and all the diagonal matrix elements = 0. 

Thus from (2.4) we can find that in the electronic eigen state i?„, the nuclei are moving in the effective potential 



Un{R) = EniR) + ^ J2 ^nkAkn , 



(2.6) 



k^n 



and transitions between the electronic states n and m are related to the non-adiabatic operator in the r.h.s. of (2.4). 
This simple observation allows us to rewrite the effective potential (2.6) as 



Un{R) = En{R) 



2to ^ 
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(2.7) 



and from this seemingly trivial expression we arrive at the following important conclusions: 

(i) corrections to En have the same order 0(7~^) as the ratio of the nuclear kinetic energy to the potential; 

(ii) off-diagonal matrix elements of the non-adiabatic perturbation operator are also small (oc 0(7^^)), and the 
fact is formulated as the so-called adiabatic theorem stating that at the adiabatic perturbations (A — > 0) there are no 
transitions between unperturbed states. 

Since the non-adiabatic effects arc characterized by the only small parameter 7"^ (the semiclassical parameter), 
the effects can be described in the frame work of semiclassical approaches (e.g. WKB or instanton ones). However, 
one has to bear in mind the main drawback problem of the Born - Oppenheimer method. Indeed the approximation 
assumes that the electronic wave functions arc real valued ones and formed the complete basis, but it is impossible 
to construct such a basis in the whole space including classically accessible and forbidden regions. 

If one relaxes the requirement to have a real valued basis, the diagonal matrix elements A„„ 7^ 0, and the effective 
adiabatic part of the Born - Oppenheimer Hamiltonian takes the form 



Un{R) 



2m 



{Wr - iAnn{R)f 



(2.8) 
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analogously to the Hamiltonian of a charged particle in a magnetic field B oc Wr x A„„. Therefore one can change 
the phases of the electronic and nuclear wave functions 

exp{ixn{R)) , ^ exp{-ixn{R)) , (2.9) 
by changing respectively the "vector potential" 

AnniR)-^ Ann{R) + ^RXn{R)- (2.10) 

Thus we confront to an important and, at times, mysterious concept of the geometrical (or Berry) phase factor that 
a quantum mechanical wave function acquires upon a cyclic evolution [38] - [47] . What is most characteristic for the 
concept of Berry phase is the existence of a continuous parameter space in which the state of the system can travel 
on a closed path. In our case the phase is determined by the non-adiabatic interaction (for more details related to the 
geometric phase for the Born - Oppenheimer systems, see e.g. the review article [48]). The phenomenom (which was 
manifested itself originally as a certain extra phase shift, appearing upon some external parameter cyclic evolution) 
has been generalized for the non-adiabatic, non-cyclic, and non- unitary cases [49], [50], although the most of the 
Berry phase applications concern the systems undergoing the adiabatic evolution (see e.g. the review article [51]). 
Note also that apart from the Berry phase some higher order corrections to the Born - Oppenheimer approximation 
(traditionally slightly confusing referred as geometric magnetism or deterministic friction, see [52]) also occur. The 
practically useful application of the Berry phase conception is the energy level displacements predicted in [53] and 
observed by NMR [54]. 

The essential physics of the phenomena can be illustrated as follows. There are two subsystems, the fast and the 
slow ones. The fast subsystem acquires the Berry phase due to the evolution of the slow subsystem. In own turn, 
there is say a feedback effect of the geometric phase on the slow subsystem. As a result the latter one is framed by a 
gauge field affecting its evolution. The gauge field produces additional forces (Lorentz-like and electric field-like ones) 
which have to be included into the classical equation of motion. In the case of stochastic external forces (e.g. from 
surrounding thermal fluctuation media), the Berry phase produces some level broadening for the fast subsystem. In 
the limit of low temperatures and strong damping, the slow subsystem dynamics can be described by the Langevin 
type equations [55]. The general message we can learn from this fact, is that the geometric phases are sources of the 
dissipative processes for LZ systems. 

Thanks to its fundamental origin, this geometric phase has attracted considerable theoretical and experimental 
attention, however, its experimentally observable consequences until now have been scarce. Therefore each opportunity 
of improving this situation is worth trying. In this respect the Born - Oppenheimer geometrical phase provides a 
unique opportunity for its observation since it must appear as a non-adiabatic contribution in a standard Bohr - 
Sommerfeld quantization rule 

S°+Xn = 2Trh, (2.11) 

where 5jJ is the adiabatic action. 

Note that care must be taken when |_E„(i?) — Em{R)\ becomes smaller with respect to the characteristic nuclear 
oscillation energy hn. It means that in the adiabatic representation (2.4) one may not consider the non-adiabatic 
interaction energy as a small perturbation. Fortunately in the limit 

\En{R) - Em{R)\ < nn 

we can start from the other limit with crossing weakly coupled diabatic states, and consider the adiabatic coupling as a 
perturbation. Thus we need the adiabatic - diabatic transformations enabling us to perform the procedure explicitely, 
which read for the wave functions 

$(i?) = exp{ieay)^{R) , (2.12) 

and for the Hamiltonians 

H = exTp(i6(jy)H exp(—i6(Ty) , (2.13) 

where (iJ, $) and {H, <&) are the adiabatic and diabatic representations respectively, (Jy is the corresponding Pauli 
matrix, and 9 is the parameter of the adiabatic - diabatic transformation (so-called adiabatic angle). 

To illustrate how it works let us consider two coupled {U12 is the coupling energy) crossing effective electronic 
potentials Ui{R) and U2{R)- The corresponding adiabatic and diabatic Hamiltonians are 

H = -^{Vnf + \{U, + U2)+ (2.14) 

Q(C/i - U2) cos2e{R) + U12 sm2e{R)^ (73 + ^ {'l^^^ ~ ^2) sin2^(i?) + U12 cos2e{R)^ cti , 
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and 

H = -|-(Vfl)2 + + U2) + J(C/i - U2)u^Ui20i , (2.15) 
im 2 2 

where 0-1,2,3 are the PauH matrices, and the adiabatic angle is chosen to ehminate the principle interaction term 
between the adiabatic states 

cot20{R) = ^^^^. (2.16) 

The adiabatic - diabatic transformation can be brought also into more elegant form [16], [56] 

{VR~iA)f^O, (2.17) 

where T is the sought for transformation matrix, and the matrix A = Ann was introduced above (see (2.5)). The 
formal solution of Eq. (2.17) can be represented as a contour integral 

f(s) = f (so)exp (- [ A{s')ds'] , (2.18) 



where sq and s stem for the initial and final points of the contour, and the solution (2.18) determines uniquely the 
transformation matrix T for the curl-less field A 

f{to) = Df{0) , (2.19) 
and the diagonal matrix D can be found from (2.17) and can be expressed in terms of the geometric phase factor 

Dkn = Skn exp(ixfc) . (2.20) 

These two relations (2.11) and (2.20) provide a complete accoimt of the non-adiabatic transitions, the cornerstone of 
the LZ problem. Besides (2.11), (2.20) show that the geometrical Born - Oppenheimer phases occur from the diabatic 
potentials crossing points, and enter quantization rules additively with the contributions from the turning points. 
Therefore, our main conclusion from this section is that non-adiabatic phenomena should be (and can be) included 
into the general semiclassical approach scheme by use of the corresponding connection matrices [57] (see also [29]) for 
the appropriate combinations of crossing and turning points of the problem. 

III. ADIABATIC PERTURBATION THEORY. 

It is almost a common student wisdom nowadays that any solution to the adiabatically time dependent Schrodinger 

equation can be represented as an expansion over the complete set of stationary (time independent) cigcn functions 
[1]. For the case under investigation (two level crossing for the electronic Hamiltonian He{r,t)), this expansion reads 

^{r,t)=c^{t)Mr)+C2{t)Mr), (3-1) 

where (pi^2 arc the stationary with respect to a nuclear motion wave functions. The time dependent Schrodinger 
equation can be rewritten exactly in the form of two first order (over time derivatives) equations for ci and C2 

where 

Hkk' = {(pk\H{t)\cPk') , k,k' = 1,2 (3.3) 

are the matrix elements for the diabatic Hamiltonian. 
The following phase transformation (see, [10], [6], [8]) 

Cfe(i) = afe(i) exp (-^ f Hkk{t)dt] (3.4) 
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reduces (3.2) to the coupled first order equations 

ihtti = Hi2a2 exp ^ J i}i2{t)dt^ , iha2 = H2iCii exp ^—i J Qi2{t)dt^ , 



wliere 



^12 = ^(-'^22 — -ffll) • 



However tlie sliglitly different piiase transformation 



Ck{t) = ^kit) exp ( ^ / {Hn + H22)dt 



keeps the second order Schrodinger like equations forms for the diabatic functions $1,2 

2 



2(i^$i / / Hii — H22 \ d 



df^ 



+ H12H21 + -^-^{Hii - H22) I <&i = , 



(3.5) 



(3.6) 



(3.7) 



(3.^ 



To clarify the mapping of this time dependent perturbation theory to the two level crossing problem, and the Born 
- Oppenheimer approach described in the previous section II, let us consider the two states Born - Oppenheimer 
equations in the diabatic representation. From (2.15) for one active space coordinate X, we have 



and 



The change of the variables 



2m rfX2 



2^ dX"^ 



(iJll -£;)$! =iJi2$2 =0, 



+ (#22 - E)^2 = #21^1 = . 



$1,2 = exp{ikoX)ci^2 , K 



2 2mE 



(3.9) 



(3.10) 



(3.11) 



transforms the two Born - Oppenheimer equations (3.9), (3.10) into the two level crossing equation (3.2) for the slow 
time dependent perturbations if one could neglect the second order derivatives {h^ /2m)d'^ci^2/dX'^ and to replace the 
time derivative by vd/dX (v = ^2Ejm is velocity). Obviously we recognize the standard semiclassical approach in 
the afore-described procedure. 

The same kind of the mapping can be performed also for the adiabatic amplitudes Ci,2(t) which are related with 
the diabatic amplitudes ci,2(i) by the adiabatic - diabatic transformation matrix depending on the adiabatic angle Q 



cos y sm t 
— sin Q cos ( 



C2(i) 



The corresponding to (3.2) set of the first order equations in the adiabatic basis 



Ci 



Hn -iO 
ie H22 



C2 



(3.12) 



(3.13) 



where the non-adiabatic coupling coefficient 6 can be related to the off-diagonal operator A12 (2.5) (or to the geomet- 
rical phase, see the previous section II) 



ie = Ax2 = i{(i>i\4>2) ■ 



(3.14) 



The transformation (3.11) enables to reduce the Born - Oppenheimer equations (for the nuclear wave functions $1,2 
in the adiabatic representation) to (3.13) if and only if: (i) in the spirit of the semiclassical approach to neglect the 
second order derivatives ; (ii) to keep only (x terms in the non-adiabatic matrix elements (i.e. neglecting higher 
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order over 1/fco contributions). Expressions (3.12) - (3.14) do allow an entry point into the adiabatic perturbation 
theory developed by Landau [1] and Dykhne [10], [11] (see also [15], [16]). We will closely follow the same method. 
One can make one step further and to find the combination of the two level system amplitudes ai,2 (3.4), (3.5) 

Y{t) = ^12^^ exp (^~^ J ^^12^*^ di + i^i2^^ exp J ni2dt^ a2 , (3.15) 



satisfying to the simple equation 



Y{t) + ^Y = 0, (3.16) 



identical to (3.8), and the both describe oscillations around the crossing point in the adiabatic potential (inverted 

adiabatic barrier). Therefore^ formally the adiabatic perturbation theory rcdTiccs the level crossing problem to the 
well known quantum mechanical phenomcnom - ovor-barrier reflection. Moreover for the latter problem the reflection 
coefficient in a full agreement with the adiabatic theorem is 1. 

Since in ID case evidently two adiabatic potentials have no reel crossing points, the crossing is possible only at 
complex values X or t 

^^12(tc) = ; Ui-U2 = ±^C/l2|f=r. . (3.17) 
In the vicinity of these points from (3.6) 

Hiaoc (^-re)^/^ (3.18) 

therefore 

^(t-Tc)3/2, (3.19) 



/ 



i.e. the crossing points are square root bifurcation points for the function ni2{t). Using Exp. (3.19) we depicted in 
the Fig. 1 the Stokes and anti-Stokes lines for the equation (3.16). The diagram shown in this figure is identical to 
that corresponding to the treated semiclassical over-barrier reflection problem with the linear turning points. The 
transition probability P12 in the main approximation is determined by the integration over the contour C{tc) going 
around the bifurcation point Tc 



P12 ^ exp (I / (J?n - H22)dt] . (3.20) 



h 



In the simplest form of the LZ problem the diabatic potentials are assumed as the linear functions of t or (what is 
the same because t = X/v) X (see Fig. 2 for the illustration) 

C/i(2) =U*±FX. (3.21) 

Putting (3.21) into the general expression for the transition probability (3.20) we find for this case 

P12 ~ exp(-27rz/) , (3.22) 



where v = Ui2/'^TT'VF is so-called Massey parameter, and the velocity v = — U'^\/m. 

Some comments about the range of the validity of the approximation seem in order. A question of primary 
importance for the LZ problem is related to the semiclassical nature of the phenomenom. To illustrate it let us note 
that for = + v'^F^X^, the Eq. (3.16) is the Weber equation with respect to the reel point X = (diabatic 
potentials crossing point). Evidently, this correspondence of two complex conjugated linear crossing points ±Tc and 
one reel X = crossing point for the Weber equation is the same as that (in the standard semiclassical treatment of 
the Schrodinger equation) for two linear and one second order turning points. Thus by the same manner as for any 
semiclassical problem, one can apply to LZ problem WKB or instanton methods. Let us compare the accuracy of the 
both approaches. If lE—lf^l » hfl (remind that f2 is a characteristic frequency of the adiabatic potentials), the WKB 
method considering for this problem two isolated linear turning points works quite well (it is the limit of koa 1, 
corresponding to the adiabatic approximation). If it is not the case we have to use the diabatic representation. 
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IV. INSTANTON METHOD IN MOMENTUM SPACE. 



This is not the place to explain the instanton method in details. However, in a stripped down version the approach 
can be viewed as follows (see [26] - [29]. [58], [59]). The recipe to find the instanton is based on the minimization of 
the functional of classical action in the space of paths connecting the minima in the upside-down potential. As it is 
well known [1] the expansion of an arbitrary wave function ^{x) in terms of the momentum eigenfunctions is simply 
a Fourier integral 



1 f+°° 



exp{ipx /h)^{p)dp . 



(4.1) 



The wave function in the momentum representation $(p) in own turn can be written down in the semiclassical form 

^p) = A{p)eM-iW{p)/n), (4.2) 
where the action W{p) is determined by a classical trajectory Xo{p) according to the definition 



dW 



(4.3) 



We use dimensionless variables: for the energy e = E/VLf), for the potential V = U/jV,o, and for the coordinate 
X = x/ao, where E and U arc corresponding dimensional values for the energy and for the potential, oq is a 
characteristic length of the problem, e.g. the tunneling distance, Qq is a characteristic frequency, e.g. the oscillation 
frequency around the potential minimum. The dimensionless momentum can be defined as 



pap 
7?l 



(4.4) 



where 7 is semiclassical parameter (remind that 7 = mfloOQ/h, where m is a mass of a particle, and we believe that 
7» !)• ^ ^ 

Introducing the semiclassical form of the wave function in momentum representation (4.2) into the standard one 
particle ID Schrodinger equation, one can transform it into the form 



2V Xo 



.l_d_ 



2 ■ 

e 

7 . 



A{P) = . 



(4.5) 



In the momentum space is a potential energy operator and it can be expanded into a semiclassical series over 1 /7 

(or what is the same, over Ti, and further wc will set everywhere Ti = 1, measuring energies in the units of frequency, 
except at some intermediate equations where it is necessary for understanding). The expansion allows us to consider 
y as a function V of two independent variables Xq and d/dP, and one can get finally 



V X, 



2dX§ 



i d \ s if dV d 

-^)=n^o) + -(^^ 



1 d'^V dXo 
+ 7: 



2dXl 



dP 



1 d^V fdXo d 



dP dP 



IcfXo 
3 dP2 



+ 



1 d'^V fdXo 



2AdX^ 



+ 



^ ,7 



d^V d^ 
dXl dP'^ 



(4.6) 



where the ellipsis represents all higher order expansion terms. 

Accordingly to the general semiclassical rules from (4.5) and (4.6) one can easily find that the first and the second 
over 7~^ order terms become identically zero, if the energy dependent trajectory Xo{P) is determined by the equation 



and so-called transport equation (TE) 



P2 + 2V{Xo) = - , 

7 



(4.7) 



dV dA 1 d^V d'^W 



dXodP 



2 dXg dP'^ 



A, 



(4.8) 
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is also satisfied. The solution of TE (4.8) can be found explicitely, and it reads 

(4.9) 



dV ^ 



dXo 

It follows from (4.9) that the scmiclassical WKB wave function (4.2) has singularities in all stationary points of the 
potential V, i.e. these points are the turning points in the momentum space. It illustrates fundamental difficulties of 
the WKB procedure, i.e. how to match the solutions which become singular on caustic lines separating manifolds in 
phase space with real and imaginary momenta. 

To illustrate also the second drawback of the WKB method let us consider the linear {V — FX) and harmonic 
(V — X'^/2) potentials. The trajectories Xo{P) can be trivially determined from (4.7). For the linear potential, 
Xo{P) is the inverted parabola with the maximum (the top of the inverted parabola) Xom — ^F/j at P = 0. The 
left and the right branches of the parabola correspond to the opposite motion directions in the classically accessible 
region Xq < Xom- The semiclassical WKB wave function in the momentum space for the linear potential 

$(P) = ^exp (-^ (eP-j^)) , (4.10) 



is the Fourier transform of the coordinate space Airy hmction. For the harmonic potential the corresponding tra- 
jectories (4.7) are the ellipses and the wave functions in the both spaces (momentum and coordinate) have the 
same functional form. It is worthwhile to note that although the WKB functions are not exact, the corresponding 
eigenvalues coincide with the exact quantum mechanical ones. 

As we have recently shown [27] - [29], one can successfully attack many important semiclassical problems by the 
instanton method. Having in mind in this section momentum space let us remind for the sake of conveniency the 
main ideas of the instanton approach. The first step of the approach derived in [58] and [59] is so-called Wick rotation 
of a phase space corresponding to a transformation to imaginary time t —it. At the transformation, the both, 
potential and kinetic, energies change their signs, and Lagrangian is replaced by Hamiltonian in the classical equation 
of motion. In the momentum space the low energy instanton wave functions can be constructed using Wick rotation 
in the momentum space (i.e. the transformation P iP) and besides the term with the energy e in (4.7) should be 
removed from this equation and taken into account in the TE (4.8). The trajectory Xq{P) in the instanton formalizm 
describes zero energy motion in the classically forbidden region of the momentum space where the wave function has 
a form 



dV 



-1/2 



HP) = ) Q{P) eM-W{P)) , (4.11) 

and additional prefactor Q{P) which has appeared in (4.11) can be represented as 

In the particular case of a linear over ID space coordinate potential {V{X) = FX) the instanton and WKB functions 

have the same form. For an arbitrary (n-order) anharmonic potential, the Schrodinger equation in the momentum 
space is reduced to the n-th order differential equation. However the derivatives of the n-th order decrease proportional 
to 7~", and therefore corresponding terms can be taken into account perturbatively. A rigorous mathematical method 
to perform this procedure (we will use in our paper) has been developed by Fedoryuk [68] - [70]. 

To illustrate the instanton approach let us consider the simplest form of the LZ problem depicted in the Fig. 3 
(see all notation in the figure caption). For the linear with arbitrary line slopes potentials in the diabatic state 
representation one has two second order coupled equations 



dX2 

2 



d^e 



= 72(a + /iX)ei =7^62; (4.13) 
= j^{a + f2X)e2=-r^vQi, 



where 9i,2 are the eigenfunctions of the corresponding states and all other notation are introduced according to the 
Fig. 3, namely 

^ = —rr- ' F = VFi F2\ , 7 = ^72— ' " = ^ — tt- ' = 2— , v = 2— . 

771(712 
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The equations (4.13) can be transformed into the momentum space (they formed a coupled linear differential equation 
system) and after that can be rewritten as one second order equation 



where we introduced 



*i = $1 exp 



.7a' 



3/2 



3 



and $1 is the Fourier transformed oi Qi, k — P/j^/a, and q{k) is a fourth order characteristic polynomial 

q{k) = X^{1 + k'^f + 2\{ik - 2u) , 

dependent of two parameters 



2(/i - /2)x/5 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



The first parameter A in the momentum representation plays the role of the new semiclassical parameter, the second 
one is the known (and already defined in (3.22)) Massey parameter. 

Fortunately all roots of the characteristic polynomial (4.16) can be found analytically quite accurately in the most 
interesting physically region of parameters. To simplify expressions (keeping complete physical content) we present 
results only for the simplest case (symmetric slopes of the diabatic potentials) /i = — /2 = /. In the classically 
forbidden region — £J > 0, a > 0, at A 1 (i.e. equivalently at a (//7)^/^), all the four roots of the polynomial 
are close to ±i 



kt 



1± 



2A 



± 



1 



2A 



— I . 



(4.18) 



In the classically accessible region (?7* — E < , a < 0), the roots are close to ±1, if A » 1 (or if —a (//7)^^^) 



kt 



1± 



4A 



1/2 



±i 



1/2 



4A 



) "'2 



-It 



4A 



1/2 



± 



1/2 



4A 



(4.19) 



(the tilde sign means that in the corresponding quantity a should be taken by its modulus, i.e. at |q|). 

Found above the set of the roots of the characteristic polynomial (4.16) in the classically forbidden (4.18) and in 
the classically accessible (4.19) regions is formally equivalent to the corresponding transition or turning points for 
the system of two potential barriers or two potential wells respectively. Thus we can use all known for these cases 
WKB and instanton results (see e.g. for all details our recent publications [29] and references herein). Since in the 
semiclassical treatment we usually concerned only with the asymptotic solutions and their connections via transition 
or turning points on the complex plane, the famous Stokes phenomenom [57], [30] of asymptotic solutions plays an 
essential role, and the distribution of the transition points (which are nothing but the zero points of the characteristic 
polynomial) and Stokes and anti-Stokes lines, dictates the phenomenom. We show all the lines emanating from linear 
turning points in Fig. 2. In the case when the roots formed a pair of close linear turning points, one can replace each 
such a pair by one second order turning points. The corresponding Stokes and anti-Stokes lines are depicted in the 
Fig. 4. 

In the classically forbidden region the instanton wave functions can be found using the roots (4.18) 



$1 = — , , . 1 exp 



Z7 



fc3 

3 



(l-ik)-'' 
*i = n — exp 



fc3 



(4.20) 



At |fc| 00 the function $J decreases proportional to \k\~^, and is reduced to the Airy function [71], [73]. In the 
vicinity of the second order turning points k = ±i, the fourth order characteristic polynomial is reduced to the second 
order one, and the equation (4.14) is reduced to the Weber equation with the known fundamental solutions [71] 



D_^{±2VX{k + i)), 
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at |fc + i| — > and 

£>_,_i(±2\/A(fc-i)), 

at |fc — i| ^0. The same procedure applied to the classically accessible region leads to the solutions 

and it is reduced to the Weber equation fundamental solutions too 

Ap(±2VIexp(i7r/4)(A; + 1)) , 

at |fc + 1| — > and 

A^-i(±2\/Aexp(i7r/4)(A; - 1)) , 

at |fc - 1| 0. 

The same solutions can be obtained for LZ problem in the two level approximation using the instanton method in 

the coordinate space. The reason for it is fairly transparent and based on the fact that for linear diabatic potentials 
the limit k ±oo corresponds to the limit x ±oo, and therefore the asymptotic of the solutions are the same in 
the momentum and in the coordinate space. 

The whole analysis can be brought into a more compact form by introducing so-called connection matrices. In 
the instanton approach we have deal with the asymptotic solutions and their connections on the complex coordinate 
plane. Thus it is important and significant to know the connection matric;es. The needed connection matrices can be 
found easily by matching of the solutions (4.20) or (4.21) at the second order turning points through the corresponding 
fundamental solutions of the Weber equation. It gives the following connection matrices 

M, = { -cos(7ri/) V2^exp(-2x)/r(i/) \ 

' \Tiiy)ex.p{2x)sm\TTu)/V2^ cos(7riv) J' ^ ' 

where X = (^ ~ (^ ~ 1/2) \ni')/2, and 

]^ ^ f -exp(-7ri>) V27rexp{-Tri))exp{-2x)/T{-w) 

\ 2r(— z;/) exp(— 7ri>/2) exp(2x)s/i(7rf>)/'\/27r exp(— 7ri>) 

where x = (-«((7r/4) + i>(l - ln(j>))) + (1/4) \ni>. 

As a note of caution at the end of this section we should also remind that initially for the linear diabatic potentials 
we had two corresponding Schrodinger equations, and each of them possesses two fundamental sohitions. Thus the 
full LZ problem is characterized by the four fundamental solutions to the left with respect to a given turning point 
asymptotic and four fundamental solutions to the right (with respect to the same turning point) asymptotic. Therefore 
generally speaking the connection matrices should be 4 x 4 ones. However owing to the symmetry of the potentials 
these 4x4 matrices have two blocks 2x2 structures for the functions $1 and <1>2- The latter ones are given in (4.22) 
and in (4.23) respectively. 




V. LZ PROBLEM FOR TWO ELECTRON STATES (INSTANTON APPROACH IN THE COORDINATE 

SPACE). 

In the previous sections II - IV we investigated LZ problem in the frame work of the adiabatic perturbation theory, 
two level approximation and momentum representation. The all three methods are equivalent and semiclassical by 
their nature. Correspondingly the approaches do work properly for the tunneling and over-barrier transmission energy 
regions, but become non-adequate within the intermediate region (of the order of 7^^^'') near the level crossing point. 
The fact is that the accuracy of these methods depends on the " renormalized" (energy dependent) semiclassical 
parameter A (4.17) which can be small in the intermediate region (A < 1 even for 7 ^ 1). To treat this region 
we have to use the coordinate space presentation, since we need to know the connection matrices for non-adiabatic 
transitions. For the latter problem the wave functions outside the level crossing point, happen to be more convenient 
(and, besides, have a more compact mathematical form) in the coordinate space. 
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A. Tunneling and over-bcirrier regions 



To move further on smoothly we start here reproducing in the coordinate space the results found in the previous 
sections for the tunneling and over-barrier energy regions. In the diabatic representation we can rewrite two second 
order LZ differential equations (4.13) as the following fourth order linear differential equation with the constant 
coefficients at the derivatives (for the sake of simplicity we consider the symmetric slope case /i = — /2 = /) 



dX 



(5.1) 



Mathematical formalism elaborated by Fedoryuk [68] - [70], we are about to recall, reduces (5.1) by a semiclassical 
substitution to a set of equations of the order 7". Characteristic polynomial for (5.1) 



F(A) = A* - 2a7'A2 - 2^^f\ + ^\a^ -v^- fX^) , 

where by the definition A = dW/dX. 

Solving the equation F{X) = perturbatively over 7"^ <C 1 we find 



Xj = A° + Uj , 



where 



and 



A° = ± 



7(a ± ^/v"^ + pX^) 



1/2 



7/ 



^3- 2 LV-j 

Four asymptotic solutions of (5.1) thus can be represented as 



{y,} = $1, $1} = {v^ + /2X2)-V4 exp 



X 



Xj{X')dX' 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



(5.6) 



and these describe a motion with imaginary momentum in the upper and lower adiabatic potentials 

2ma^ 



(U^ -E)= 72(a ± VV+72X2) . 



The subscripts in (5.6) corresponds to the upper or lower adiabatic levels, and the superscripts are referred to the 
sign of the action. 

Before turning on the connection matrices let us use the substitution 



$1 = exp{KX)(j), 
and choose n value to vanish the first derivative in (5.1), i.e. 

- -f'^an - ^7^ = . 

At a > 3(//47)^/^ one can expand the roots of (5.8) in terms of the parameter 

47 

Thus we find 



3\/3 ■ 



(5.7) 



(5.8) 



(5.9) 



(5.10) 
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At the condition (5.8) the coefRcients at the fourth and at the third order derivatives in (5.1) are small (proportional 
to S and to VS respectively) and the fourth order equation (5.1) can be rewritten as two second order Weber equations 
with the solutions 

where 

,._.4_.(,_-)„..|_„(.4),,.,.(i^)'"(,4). 

The leading terms of these solutions are the same as found in the previous section IV. But the Fedoryuk method we 
used, gives us also in the tunneling region (5.8) the higher order over the parameter i5 corrections. 

In the over-barrier energy region where a < — 3(//47)^/^, the roots of the equation (5.8) are complex - conjugated 



and S plays the role of the small parameter in this region 




/ 



47|a|3/2 • 



(5.12) 



(5.13) 



Again as above for the tunneling region, the coefficients at the higher order derivatives are small, and therefore, the 
function (j) (5.7) satisfies the Weber equation with the fundamental solutions 

-Dp(i,2) (/3(i_2)^) ) 

where 

= + +w + ^ \ ,f=i^+w[l-^\ , (5.14) 




' 2 j-2 \ 1/4 / 2 f2 \ 1/4 

Pi = exp(z7r/4) ( ) , p2 = exp(-i37r/4) ' ^ ^ 



a\ j \ |q; 

Like it was for the tunneling region (5.11), the leading terms of the expansion (5.14) coincide with the results found 
in the previous sections, but from (5.14) we are able to compute the corrections to the main terms. 

Now we are in the position to find the connection matrices. To do it for the tunneling region we have to establish 
the correspondence between the solutions of the fourth order differential equation (5.1) and those for the localized 
in the left (L) and in the right (i?) wells states. In the case a » for the diabatic potentials, the action can be 
computed starting from the both wells {R and V) 

^W^ ~ 7Wo^ + koX + '—X^ , jW'^ ~ jW^ - koX + ^—X^ , (5.15) 

where ko = {2ma'^{U* - E)lh^fl'^ = 7\/a is imaginary momentum, and W^q ' are the actions computed from an 
arbitrary distant point in the L o r in the R wel ls, respectively to the point X = 0. From the other hand in the 
adiabatic potentials = U"^ ± a/ [7^2 + a^pX'^ the corresponding actions can be represented 

^W^ -^W^ = koX±^X'^sign{X). (5.16) 

Explicitely comparing the semiclassical wave functions in the both representations (adiabatic and diabatic ones) it is 

easy to see that the adiabatic functions in the potential /7~ coincide with the diabatic functions for localized L and 
R states at X < and X > respectively. The adiabatic functions for the upper potential correspond to the 
tails of the diabatic wave functions localized in the opposite wells. Therefore in the level crossing region the L/R 
diabatic functions are transformed into the R/L functions, and the interaction entangles the diabatic states with the 
same sign of k^X. Thus we have only four non-zero amplitudes of the following transitions 
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Recalling that 



1/2 «2 

~ koX ± —X^ ± -(1 -\niy). 



(5.16 



we come to the conchision that the quantum sohitions (5.11), vaUd in the vicinity of the level crossing point asymp- 
totically, match smoothly increasing and decreasing solutions (5.6), and it leads to the Landau description [1] of the 
level crossing transitions depicted in the Fig. 5. 

Using expressions (4.22), (4.23) relating the fundamental solutions of the Weber equation, we can find the corre- 
sponding to (5.17) 4x4 connection matrix 



(5.19) 



27rexp(-2x)/r(j/) 



cos(7ri/) 



where as above 





T{v) exp(2x) sm^{ni') 

cos(7ri^) 




X = 





■ cos(7rz/) 

27rexp(-2x)/r(i^) 




— cos(7ri/) 



r(i/) exp(2x) sin2(7ri/)/V27r 



The matrix (5.19) has a 2 x 2 block structure, and each of the identical blocks connects increasing and decreasing 
diabatic solutions. However these diagonal blocks do not correspond to the L — R transitions separately for the lower 
and upper adiabatic potentials. Indeed the corresponding to these transitions 2x2 matrix is 



2ir exp(— 2x) /r(z^) — cos(7ri/) 

cos(7rz^) r(i/) exp(2x) sin^(7r^')/\/27r 



(5.20) 



In the diabatic limit (i.e. v ^ 0) the diagonal matrix elements are small (oc i/^/^, and v^^"^ respectively), and the 
off-diagonal elements approach to ±1, as it should be since by the definition there are no transitions between the 
diabatic potentials. 

In the adiabatic limit v ^ 1, the diagonal matrix elements tend to 1, and it means that the decreasing L solution 
transits only into the increasing R solution, and vice versa. Thus the connection matrix in the tunneling region 
depends only on the Massey parameter v. One has to bear in mind here that the blocks of the 4x4 connection matrix 
(5.19) correspond to the two isolated second order turning points with the Stokes constant (see e.g. [29]) 



To 



v^27r 



exp(-2x) . 



(5.21) 



Analogously one can study the over-barrier region. Repeating again the procedure described above for the tunneling 
region (with evident replacements ko —iko and 1/3"^) we end up with the following 4x4 connection matrix 



U 



27rexp(-2x)/r(-iz/) 

2r(— iz/) exp(— TTi/) exp(2x) sinh(7ri/)/-\/27r 

exp(— TTi/) 

exp(— TTI/) 

— exp(— TTi^) 

— exp(— TTJ^) 

27rexp(-2x)/r(-iz/) 

2r(— iz/) exp(2x) exp(— ttz^) sinh(7rz/) / \/2tt 



(5.22) 



where 



(J + Kl-lni^)) + J(7rz^ + ln^). 



I /TT 

2 



(5.23) 
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The same manner (as it was already mentioned for the tunneHng region), the blocks in (5.22) correspond to the two 
isolated second order turning points with the Stokes constant [29] 

^2 = f^e^P(-2x). (5.24) 

Thus we arrive at the important conclusion that the main peculiarity of the LZ level crossing (in comparison with 
the standard say one-potential problems) is that the second order turning points characterizing the diabatic levels 
crossing for the LZ problem, possesses the different Stokes constants T2 (5.21) and T2 (5.24) in the tunneling and in 
the over-barrier regions. 



B. Intermediate energy region 



We can now reap the fruits of the previous subsection consideration efforts. First, let us note that from the relations 
(5.11) and (5.14) one can see that when the energy approaches to the top of the barrier, the exponents p^^^ and p^*^ 
of the parabolic cylinder functions are increased and thus, more and more deviated from the value prescribed by the 
Massey parameter z/. Second, increasing of /3(j) upon \a\ decreasing, decreases the values of \X\ where the asymptotic 
smooth matching of the solutions should be performed. For 5^0 these \X\ values are located deeply in the classically 
forbidden region, where the potentials are close to the diabatic potentials, while for S > 2-\/3/3, these coordinates \X\ 
are of the order of the quantum zero-point oscillation amplitudes, and therefore to find the solution in this region, we 
have to use the adiabatic representation. 

These two simple observations give us a conjecture how to treat LZ problem in the intermediate energy region. To do 
it first of all we should find the energetical "window" for the intermediate region. It is convenient to chose the adiabatic 
potential frequency O = F/^/mU\2 as the energy scale, and in terms of this scale the inequality \a\ < 3|//(47)|^/^ 
reads as 



By the other words the characteristic interaction energy at the intermediate region boundaries does not depend on 
Ui2- However, the positions of the linear turning points \X*\ corresponding to the energies U* ± depend on the 
ratio Ui2/Ui2- These points are located inside or outside of the interval [— ao7~^/^ +ao7~"^/^] at U12/UI2 < 1 and at 
U12/U12 > 1, respectively. Accordingly for the both cases the matching conditions in the intermediate energy region 
are different. In the former case for the asymptotic matching region the potentials can be reasonably approximated by 
parabola, and therc^fore we should work with the Weber equations, and for the latter case the matching are performed 
in the region where the potentials are linear ones, thus the equations are reduced to the Airy equations. 

Let us investigate first the case U12/U12 > 1- Using Born-Oppenheimer approach of the section II, the Schrodinger 
equations in the adiabatic representation with the accuracy up to 7"^ are decoupled for the wave functions ^± 



dX^ 



+ 7^(0! ± ^v^ + /2X2)5'± = . (5.26) 



The equations (5.26) at \X\ < v/f are reduced to the Weber equations with the fundamental solutions 

D_i/2-q,{±V^X), and D_i/2+iq^{±exp{-iTT/4:)y/2jX), where 

V + a V — a ,_ 

91 = , Q2 = , (5-27) 

do not depend on the Massey parameter v. Two real solutions of (5.26) correspond to the upper adiabatic potential 
(classically forbidden region), and two complex solutions correspond to the classically allowed motion under the lower 

adiabatic potential. 

The argument of the Weber functions is oc X^/^, and at the condition X < v/f their asymptotic expansions 
determine the interval where the matching should be done 

7-^/^(^)'''>7-^/^ (5.28) 



Thus this inequality (5.28) can be fulfilled only at U12/UI2 > Ij when the intermediate region is sufficiently broad in 
comparison with fl. In this case the exponents qi , (72 (5.27) are large, and our aim is to find explicitely the asymptotic 
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expansions of the solutions in this case. Wc will closely follow the method we borrowed from Olver paper [72] (see 
also his monograph [73]), which is in fact an expansion over small parameters l/|gi| (where are the exponents 
(5.27)) of the fundamental Weber solutions, and it leads to the following asymptotic solution to Eq. (5.26) at X > 



1/2, 



(5.29) 



where Y± = ^/v zL a + X"^. Using the known relation between the fundamental solutions of the Weber equation [71], 
[73] 



^/i(^) = exp(-i7r/u)£'^(z) 



2tt 



r(-M) 



exp 



-ITT- 



D. 



{iz), 



we can find two other (complimentary to (5.29) solutions 

*i(X) = y-^/' ( - sin(7rgi)(X + F+)-«^ exp(-7Xy+) + exp(-2xi) 



27r 



r((l/2) + gi) 



{X + Y+reMl^y+) (5-30) 



and 



^liX) 



(5.31) 



1/2 



iexp{-Trq2){X + Y-Y"^ exp{ijXY-) + exp(-2x2) 



27r 



r((l/2) - iq2) 



where we introduce the notation 



Xi 



X2 



iq2 



2 



(X + y-)*«= exp(-nxy_) 



In 52 



Not surprisingly but it is worth to noting that these solutions (5.29) - (5.31) can be represented as a linear combination 
of the semiclassical solutions (5.6) with the coefficients 



cos 2^(1^2) 



X 



Vv±a + X^ 



(5.32) 



These energy dependent angles ^(1,2) coincide with the adiabatic angles (see (2.12) and (3.12)) introduced above in 
sections II and III in the level crossing point at a = 0, and < v, and the both angles aquire only slightly different 
values over the whole intermediate region jaj < v. 

Now we can find all needed connection matrices for these functions. Although the calculation is straitforward it 
deserves some precaution (e.g. the X-dependent matrices have different functional form at the positive and negative 
X). At X > we get 



cos 62 
-i exp(— 7r(72) cos 6*2 







2^exp(-2x2) cos6'2/r((l/2) - 152) 





(5.33) 






sin6'i 
sin(7rQ'i) sin 01 








27r exp(-2xi) sin0i/r((l/2) + q^) 



and for X < it reads as 



27rexp(-2x2) cos6'2/r((l/2) - iq2) -i exp{-'rrq2) cos 82 
cos 62 






(5.34) 






sin6'iA/27rexp(-2xi)/r((l/2) + 91) - sin(7r(7i) sin 6*1 
sin 6*1 
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The product of the inverse to (5.33) matrix and the matrix (5.34) determines the connection matrix, we have sought for, 
to relate the semiclassical solutions in the intermediate energy region (cf. (5.20) and (5.22) presenting the connection 
matrices for the tunneling and over-barrier energy regions). Performing this simple algebra one ends up with 



Ucross — 



\/2^exp(-2x2)/r((l/2) - iq2) texp(-7rg2) 

-i exp{-nq2) 2 exp(2x2)r((l/2) - iq2) cosh(7rg2) /r qr^ 








V2^exp(-2xi)/r((l/2) + gi sinTrgi 

-sinTT^i cos^(7rg'i)r((l/2) -|-g'i)exp(2xi) 

The matrix (5.35) has two 2x2 blocks structure, the same as the connection matrices (5.19) and (5.22) for the 
tunneling and over-barrier regions. However, unlike (5.19), (5.22) describing the transitions between the diabatic 
states, the matrix (5.35) corresponds to the transitions between the adiabatic states. Indeed, at a strong level 
coupling {Ui2 > U12) the eigenf unctions are close to the adiabatic functions and only non-adiabatic perturbations 
induce the transitions. Respectively, the off-diagonal matrix elements in (5.35), having meaning of the probability 
to keep the same diabatic state after the transition, are zero. The block with the real - valued matrix elements 
corresponds to the minimum of the upper adiabatic potential, i.e. it is to the isolated second order turning point, 
where [29] 

1 U*-E + 1/12 ^ 1 

'^^ + 2 = n + 2- (^-^^^ 

The complex - valued block is associated with the maximum of the lower adiabatic potential, and analogously to 
(5.36) one can find the following relation for the turning point 

In the case of weak level coupling, namely at \U* — E\ < UI2 and Ux2 < U12, for the intermediate energy region, the 
adiabatic potentials everywhere (except a small neighbourhood \X\ < v/f — > of the level crossing point) can be 
linearized, i.e. represented as a ± f\X\, and the asymptotic solutions (5.6) are reduced to a linear combination of the 
following functions 

$± oc (/|X|)-V2 eMH+signX) , oc (/|X|)-V2 eM^^-signX) , ^± = §j{f\X\ ± af^^ . (5.38) 

Now all needed matrix elements can be calculated in the frame work of the Landau perturbation theory [1], which 
can be formulated to avoid divergency of the parameter at a — > in terms of the dimensionless variables 

a = 3.2-^/3C{^^^3.2-4/3t^_ 
U12 U^2 

The results of our analysis is shown in Fig. 6. The tunneling and the over-barrier regions are separated from the 
intermediate energy region by the lines |C^i2~^l = ^^12- own turn the intermediate region is also separated into two 
parts by the line v = v* = 0.325, where ly* is the value of the Massoy parameter ly at U12/UI2 = 1 and \U* — E\ = [/j*2- 
In the V < V* region the perturbation theory is the adequate tool for the problem, and the transition matrix elements 
are proportional to ?7i2/C^i2- At v > v* one can use the connection matrix (5.35). To illustrate the accuracy of the 
approxiniMtions wc have computed the matrix element Afn. The results are shown in the Fig. 7. Our computations 
demonstrate quite good precision, secured up to two stable digits. The accuracy of the results on the boundaries 
between the intermediate and over-barrier or tunneling regions is not worse than 3 — 5%, and can be even easily 
improved using interpolation approaches. 



VI. SCATTERING MATRIX. 

LZ kind of phenomena can be considered as (and applied to) scattering processes. Found in the previous section 
V expressions for 4 x 4 connection matrices can be used to calculate the scattering operator (or matrix) S, which 
converts an ingoing wave into an outgoing one. 
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Let us consider first the over-barrier region for two linear potentials crossing problem. For this case besides 
the crossing point which we chose as X = 0, there are two linear (first order) turning points Xq = ±\a\/f (each 
turning point for each of the diabatic potentials, we designate by L and R). The scattering matrix which relates the 
asymptotic solutions at X <C —Xq and at X Xq is the product of the 4x4 connection matrix (5.22) calculated in 
the section V (subsection VA) and known [57] (see also [29]) two semiclassical connection matrices describing wave 
function evolution from the turning point —Xq to the crossing point 0, and from this point to the turning point +Xo, 
respectively. Thus we end up with 2x2 matrix having the following block matrix elements 



Tn 



if 



expi{^—(po) 

exp — (po) 



(6.1) 



= t;, = (1 - A^f) eMhW*/2) 

T22 = Aif 



i -1/2 
■exp—i^W* (i/2)exp— 17!^* 



2 cos{-fW* -{(j>^ - sin(7T^* - (0 - 0o)) 

sin(7W^* -((/.- 0o)) (1/4) cos(7W^* -(</.- (/.q)) 



X (see 



where Aif = (1 — cxp(— tti^))"'^/^ is the LZ amplitude of the transition between the diabatic states, (j) — 4 
(5.23)), and W* is the action between the linear turning points. 

The diagonal elements of (6.1), proportional to the transition amplitude A^f, describe the propagating waves (i.e. 
the solutions of the Schroedinger equation in the lower adiabatic potential), and oscillating blocks correspond to the 
solutions in the upper adiabatic potential. Off-diagonal blocks proportional to the probability to keep unchanged the 
initial diabatic states, describe the waves reflected from the linear turning points. Interesting for physical applications 
reflection R and transmission T coefficients can be found from (6.1) by a straightforward calculation, and the results 
read as 

R = - Alf)[Al- exp(*7T^* - 2i{c^ - 0o)) + exp(-z7T4^*)]-i ; (6.2) 
T = 2A,f cos(7iy* - (0 - </)o))[A2^ exp(z7Ty* - 2i(0 - 0o)) + exp(-z7M^*)]-i . 

The poles of the scattering matrix can be also easily found from (6.1), and the corresponding resonance condition is 



cos(2(7l^* -{(}>- <po))) = ~(^^-\ exp(-27rz/)^ (1 - exp(-27rz/))-^/^ . 
The action in the resonance points is complex- valued 

Rei'yW* -{(f>- <^o)) =(n+^^7r; Im{'yW* -{(f>- M) = ln(l - exp(-27ri/)) . 



(6.3) 



(6.4) 



The poles of the scattering matrix are placed on the lower complex E half plane at the vertical lines corresponding 
to conventional Bohr - Sommerfeld quantization rules (7W* = 7r(n -I- (1/2))) for the upper adiabatic potential. In 
the diabatic limit {v — > 0) the imaginary part of the pole positions tends to infinity, and in the adiabatic limit 
{1/ 00) the poles move to the real axis. Thus we see that the eigenstates of the upper adiabatic potential are always 
quasistationary ones. The resonance widths are determined by the residues of the scattering matrix elements at the 
poles, and it can be shown the resonance widths are monotonically decreasing functions of ly. In the Fig. 8 we show 
the energy dependence of the transmission coefficient T. In the diabatic limit T — > 0, and it is increased when U12 is 
increased, and in the over-barrier region there appear the resonances with the widths r„ increasing with the energy 
increase, since in these conditions the Masscy parameter is decreased and r„ oc exp(— 27rz/). 

We illustrate the energy dependence of the transmitted wave phase in the Fig. 9. In accord to the general scattering 
theory [1] , there are tt- jumps of the phase at each quasi-discrete energy level of the upper adiabatic potential. At 
U12/U12 < 1, the resonance widths are of the order of the inter- level spacings. The amplitudes of the decaying 
solutions (localized in the well formed by the upper adiabatic potential) are increased near the resonances, and this 
behavior is depicted in the Fig. 10. One note of primary importance concerning the issue is that the information 
about decaying solutions existing in the 4x4 connection matrix, e.g. (5.22) is lost when we use 2x2 scattering 
matrix (6.1). 

Except for a slight natural modification the presented above results, one can find the scattering matrix for the 
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tunneling region merely by recapitulating the already derived expressions. Thus instead of the matrix (6.1) we get 



Til = 



(l/4)Miiexp(--7T4^*) + M22Cxp(7l'F*) i((l/4)Mii cxp(-7T4^*) - M22 exp(7VK*)) 
-i((l/4)Mii exp(-7lV*) - M22 cxp(7VK*)) (l/4)Mii exp{-jW*) + M22 expi-jW*) 

l/2exp{--fW*) 



(6.5) 



T12 = T2*i = cos(7rz/) exp(i7VF72) 
T22 



-1 (i/2)exp(-7W^*) 



Mil 
M22 



where Mn and M22 are the corresponding matrix elements from (5.19). 
The same as above we compute the reflection and transmission coefficients 



1 



1 



R=-i{ exp(7T4^*) - -Mf^ exp(-7T4^*) exp{jW*) + -M(^ exp(-7VF*) 



T = Mil ( exp(7W^*) + -Mfi exp(-7W^*) 



(6.6) 



In the intermediate energy region the only block matrix element Tn requires the special calculations taking into 
account the contributions from the complex turning points 



^11 



27rexp(-7r(j2/2)/r((l/2) - iq2) i cxp(-7rg2) 

-i exp(-7rg2) 2r((l/2) - iq2) cxp{-TTq2/2)ch{Tiq2)lV2^ 



(6.7) 



and all other matrix elements are the same as in the matrix (5.34). Finally we find also in the intermediate energy 
region the reflection and the transmission coeflicients 



R = 



exp(-7rg2) 



-y/l + exp(-27rg2^ 
where (f) = <ffl'[r((l/2) — 152)]. 



:expH(<^-(7r/2))] ;r = 



^Jl + exp{-2'Kq2) 



exp(— 



(6.^ 



VII. QUANTIZATION RULES FOR CROSSING DIABATIC POTENTIALS 

In spite of the fact that instanton trajectories are rather simple objects, and can be relatively easy found analytically, 
calculations of the quantization rules within the instanton approach are rather involved and intricate and require the 
knowledge of the scattering matrix and all connection matrices, we have calculated in the previous sections. In this 
section we apply this machinery to find the quantization rules for the crossing diabatic potentials depicted in Fig. 
11. Depending on the Massey parameter the shown on the figure situations exhaust all cases practically relevant for 
spectroscopy of non-rigid molecules (symmetric or asymmetric double - well and decay potentials). 

Within the instanton approach the quantization rule can be formulated as a condition that the amplitudes of 
exponentially increasing at X > 0, and X < Q, respectively, solutions <1>^, must be vanished. Taking into account 
that W£ = (the actions in the corresponding wells of the lower adiabatic potential) and using found above in the 
section IV the connection matrix (5.19) , the quantization rule for this case is 

tan(7iy2)=±^exp(7iy^), (7.1) 

where VF^ is the action in the barrier formed in the lower adiabatic potential, and p = Uu is the corresponding matrix 
element of the connection matrix (5.19). 

Only the factor 1/p varying from to 1 in the diabatic and in the adiabatic limits, respectively, makes this 

quantization condition (7.1) different from the well known [1] quantization rule for the symmetric double- well potential. 
Correspondingly, the tunneling splitting at finite values of the Massey parameter i/ can be represented as a product 

A„ = AOp(i^), (7.2) 

of the tunneling splitting A° in the adiabatic potential and the factor 

Pi'^) = f^.l''-^'^'^ ^M-'^) , (7.3) 



20 



associated with the transition amplitudes between the diabatic potentials in the crossing region. 

It is particularly instructive to consider (7.1) as the standard [1] Bohr-Sommerfeld quantization rule, where in the 
r.h.s. the both, geometrical Lpn and the tunneling x„ phases are included additively. In the adiabatic limit when 
p{y) 1, we find that (pn ^ and (7.1) is reduced to the quantization of the symmetric double-well potential. In 
the diabatic limit (fin = —Xn and the geometric phase compensates the tunneling one. The physical argument leading 
to this compensation may be easily rationalized as follows. Indeed, at the reflection in the crossing point X = 0, the 
trajectories in the classically forbidden energy region are the same as those for the tunneling region but with a phase 
shift TT. 

We focus now on the quantization rules for the over-barrier energy region. Closely following the consideration 
performed above for the tunneling region (replacing the connection matrix (5.19) by the matrix (5.22), and making 
some other self-evident replacements) we end up after some tedious algebra with the quantization rule 

(1 - exp(-27r!/)) cos(27PF£ + ((/>- ^o)) cos{-fW* -{(t>- </>o)) + exp(-27ri^) cos^ {^iWl + ) = , (7.4) 



where W* is the action in the well formed by the upper adiabatic potential, and ^ — 0o = X is determined according 
to (5.23). Prom the Eq. (7.4) follows that the eigenstates are determined by the parameter 

B = ^"P(''":) , . (7.5) 
1 - exp(-27rj^) ^ ' 

In the diabatic limit v ^ {), and therefore, B l/{2'Ky) in (7.4) the main contribution is due to the second term, 
and it leads to a splitting of degenerate levels in the diabatic potentials. Moreover since 



W*\ ( 1 
7( W^L + ^ 1 =7r( n+-±2i/sin 



(7.6) 



the splitting increases when the Massey parameter v increases, and it is an oscillating function of the interaction Z7i2. 

In the adiabatic limit, when v oo, (/) — 0o ^ 0, and, therefore, from (7.5) B ~ exp(— 27ri^), the main contribution 
to (7.4) comes from the first term which determines the quantization rule for the upper one-well potential and for 
the lower double-well potential in the over-barrier energy region, and in this limit the parameter B plays a role of 
the tunneling transition matrix element. For B smaller than nearest level spacings for the lower and for the upper 
potentials, one can find from (7.4) two sets of quantization rules leading to two sets of independent energy levels 

iW* = TT (^m + ; I^Wl = TT (^712 + . (7.7) 

Since the eigenstatc energy level displacements depend on XJyi the resonances can occur at certain values of this 
parameter, where the independent quantization rules (7.7) are not correct any more. The widths of these resonances 
are proportional to exp(— 27r!/) and therefore are strongly diminished upon the Massey parameter v increase. This 
behavior is easily understood, since in the limit the wave functions of the excited states for the lower potential are 
delocalized, and their amplitudes in the localization regions for the low-energy states of the upper potential, are very 
small. 

More tricky task is to derive the quantization rule in the intermediate energy region. One has to use the connection 
matrix (5.35), and to bear in mind the contributions from the imaginary turning points. Nevertheless, finally the 
quantization rule can be written in the simple and compact form as 

cos(27T4^2) = - exp(-7rg2) , (7.8) 

where = 7(1; — a)/2 is determined by the relation (5.27). 

It is useful to illustrate the essence of the given above general result by simples (but yet non trivial) examples. 
First, let us consider two identical parabolic potentials with their minima at X = ±1 and with the coupling which 
does not depend on X. Since the symmetry, the solutions of the Hamiltonian can be represented as symmetric and 
antisymmetric combinations of the localized functions 

= ^($L ± ^r) . (7.9) 

The functions are orthogonal, and, besides, two sets of the functions (^f^ , ^'q )> and (^'J , ^^) (where the subscripts 
and e stand for the ground and for the first excited states respectively) correspond to the two possible kinds of level 
crossings. 
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In Fig. 12 we depicted schematically the dependence of the level positions on the coupling Ui2. In the energy region 
E < U* + Ui2 where only there exist the discrete levels of the lower adiabatic potentials, there are the pairs of the 
alternating parity levels (^'+ , ), and (^'^J" , ^'j). The tunneling splittings are increased monotonically since the 
Massey parameter v is increased, and the barrier is decreased with Ui2. The same level and parity classification is 
remained correct for the energy region above the barrier of the lower adiabatic potential where the spectrum becomes 
almost equidistant one. However, in the over-barrier region, the resonances are occurred between the levels of the 
same parity, and this sequence of the odd and of the even levels is broken, and level displacements are not monotonic 
functions of Ui2- Some of the levels of different parities can be mutually crossed. For the upper adiabatic potential the 
level sequence is opposite to this for the lower adiabatic potential. Note that wc checked the results of our semiclassical 
approach and found remarkably good agreement with the numerical quantum diagonalization. 

The second instructive example treats the one-well and linear diabatic potentials crossing. It leads to the lower 
adiabatic decay potential and to the upper one- well adiabatic potential. The quantization rules in this case correspond 
to the vanishing amplitudes for the exponentially increasing solutions when X ^ —oo, and besides one has to require 
that no waves propagating from the region of infinite motion, i.e. a.t X > 1/2. Performing the same as above procedure 
we find that in the tunneling energy region, the eigenstates are the roots of the following equation 

tan(7W*) = -i-i-y exp(27W^^) , (7.10) 

with the same as above notation. 

To proceed further it is convenient to introduce the complex action to describe the quasi - stationary states 

^Wl=n(^-i^], (7.11) 



where evidently fl = dW^/dE does depend on E. From (7.11) the reel and imaginary parts of the quantized 
eigenstates are 

En=^(n+^ ■,Tn= p2(i.)il exp(-27l^|) . (7.12) 

This relation (7.12) describes the non-adiabatic tunneling decay of the quasi-stationary states of the lower adiabatic 
potential. The same as we already got for the two parabolic potentials crossing (7.2), here the tunneling and the 
adiabatic factors are entering decay rate multiplicatively. Since the decay rate is proportional to the square of the 
tunnehng matrix element, r„ oc p^{i') as it should be. 
In the over-barrier energy region the quantization rule is 

(1 - exp(-27rz/) exp[-i(7l^2 + + cos(7H^* - (/-o) + (7.13) 
exp(-27rz/) eyL^{-i'yW* /2) cos i^Wl + ) = , 



and the actions depend on the energy E as 



U* + Ui2 E_ 



(7.14) 



where Vl and Q.\ are -E-dependent frequencies of the diabatic and the upper adiabatic potentials. 
In the diabatic limit the decay rate is proportional to the Massey parameter and has a form 



r„~7rz/cos2(7VF-(^-h(/)o), (7.15) 

and in the opposite, adiabatic, limit the decay rate is 

r„ ~ exp(-27r!.)(l - sm{2-iWl + <^ - <^o)) • (7.16) 

In the both limits the decay rate is the oscillating function of Ui2. Wc illustrate the dependence T{Ui2) for the 
crossing diabatic potentials U\ = [l + X)"^ /2 and U2 = (1/2) — X in the Fig. 13. Note that while the tunneling decay 
rate of the low-energy states is increased monotonically with the Massey parameter v, the decay rate of the highly 
excited states goes to zero in the both (diabatic and adiabatic) limits. Besides there are certain characteristic values 
of Ui2 when the r.h.s. of (7.15) or (7.16) equal to zero and therefore r„ = 0. 
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Last and more general example we consider in this section, describes two non-symmetric potentials crossing at 
X = point: 

Ui = ^(1 + Xf ; U2 = ^{X^ - 2bX + b) . (7.17) 

In a certain sense it is the generic case, and when the parameter b entering the potential (7.17) is varied from 1 to oo, 
wc recover the two particular examples considered above, and come from two identical parabolic potentials to the case 
one-well and linear diabatic potentials crossing. This kind of the potential U2 was investigated recently by two of the 
authors (V.B. and E.K) [64] aiming to study crossover behavior from coherent to incoherent tunneling upon increase 
of the parameter 6, the larger is this parameter b, the larger will be the density of final states. The criterion for 
coherent-incoherent crossover behavior found in [64] based on comparison of the transition matrix elements and the 
inter level spacings in the final state. The analogous criterion should hold for LZ level crossing problem, however in 
the latter case the tunneling transition matrix elements has to be multiplied by the small adiabatic factor. Therefore 
the coherent - incoherent tunneling crossover region moves to the more dense density of final states, and the larger 
U12 is the smaller will be the region for incoherent timneling. 

Quite different situation occurs for highly excited states. In the diabatic limit, the transition matrix element is 
increased with the Massey parameter i/, and therefore at a given b value, the system moves to more incoherent 
behavior. In the adiabatic limit, the transition matrix element is exponentially small, and coherence of the inter-well 
transitions should be restored. However, since the matrix elements are oscillating functions of U12 for the intermediate 
range of this coupling {U12) coherent - incoherent tunneling rates are also non-monotonically varying functions. To 
illustrate these unusual phenomena we show in Fig. 14 time dependence of the survival probability P for the initially 
prepared localized in the left well state n = 0. 

VIII. CONCLUSION. 

In this paper we have challenged again the very basic subject - LZ problem. Currently there are about 100 
publications per year related to LZ problem. Clearly it is impossible to give a complete analysis of what is being 
achieved in this field. Our aim, therefore, was only to show some recent trends and our new results, to help newcomers 
and specialists in finding cross-references between the many physical phenomena related to LZ problem. The problem 
was first addresses a long ago, and many already classical results are known now from the textbooks [1] , [37] . Although 
exact quantum- mechanical calculations are still prohibitively difficult, many important results have been obtained in 
the frame work of the WKB approach [1] -[65]. The accuracy of the modified WKB methods can be improved 
considerably, note for example [30] where the authors have included into the standard WKB method additionally 
a special type of trajectories on the complex phase plane, along which the semiclassical motion is described by the 
Weber functions. This method ascending to Landau [1] is equivalent to the appropriate choice of the integration 
path around the turning point, and it appears to be quite accurate for the tunneling and over-barrier regions, where 
the characteristic fourth order polynomial (see (4.16)) can be reduced to the second order one (two pairs of roots 
are nearly degenerated). However, even in this case there are some non- negligible corrections found in the papers 
[23] - [25]. In the intermediate energy region, where all 4 roots are noticeably different, the method becomes invalid. 
Besides the choice of these additional special trajectories (which one has to include to improve the accuracy of the 
WKB method near the barrier top) depends on the detail form of the potential far from the top, and therefore for 
each particular case the non-universal procedure should be perform from the very beginning. 

We believe we are the first to explicitely addresses the question on the behavior in the intermediate energy region. 
In all previous publications this region was considered as a very narrow and insignificant one, or in the best case the 
results were obtained by a simple interpolation from the tunneling (with monotonic decay of the transition probability) 
to the over-barrier (with oscillating behavior) regions. The fact is that classical trajectories can be separated into two 
classes: "localized" and " delocalized" in the following sense. If energy is close enough to the minimum or maximum 
of the potentials, the trajectories could be called confined, since they are determined by the universal features of the 
potentials in the vicinity of these extremal points. Evidently it is not the case for the intermediate energy region. 
In the present study we have found that, contrary to the common belief, the instanton trajectory is a rather simple 
object and can be explicitely computed even for the intermediate energy region. 

Within the framework of the instanton approach we present a full and unified description of ID LZ problem, which 
is very often can be quite reasonable approximation for real systems. Because different approaches have been proposed 
to study the LZ problem we develop an uniform and systematic procedure for handling the problem. We reproduced all 
known results for tunneling and over-barrier regions, and studied as well the intermediate energy region. Specifically 
we applied our approach to the Born - Oppenheimer scheme, formulated the instanton method in the momentum 
space, and presented all details of the LZ problem for two electronic states using also the instanton description of 
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the LZ problem in the coordinate space. Neglecting higher order space derivatives we find asymptotic solutions, and 
using adiabatic - diabatic transformation we match the solutions in the intermediate region. Based on these results 
we derived the complete scattering matrix for the LZ problem, the quantization rules for crossing diabatic potentials. 
Our results can be applied to several models of level crossings which are relevant for the interpretation and description 
of experimental data on spectroscopy of non-rigid molecules and on other systems undergoing crossing and relaxation 
phenomena. 

Note also, that in spite of the fairly long history of the LZ phenomena, the study is still in an accelerating stage, 
and a number of questions remain to be clarified (let us mentioned only few new features of the phenomena attracted 
attention recently, like LZ interferometry for qubits [74] , LZ theory for Bose - Einstein condensates [75] , multi-particle 
and multi- level LZ problems [76], [77], [78], [79]). Much of the excitement arises from the possibility of discovering 
novel physics beyond say the semiclassical paradigms discussed here. For example, we found in the sections II and III 
that the wave functions of nuclei moving along the periodic orbits acquire geometrical phases (the effect is analogous 
to the Aharonov - Bohm effect [38] , but in our case it has nothing to do with external magnetic fields and is related to 
the non- adiabatic interactions). The relation between the both phenomena (the geometrical phases and the periodic 
orbits) can be established using Lagrangian (instead of Hamiltonian) formulation of the problem, which enables to 
take into account explicitely, using propagator technique, [34], [35], [36], time dependence of the adiabatic process 
under consideration (see also, e.g., [4], [43]). However, a proper handling of these aspects is beyond the scope of 
our work. Further experimental and theoretical investigations are required for revealing the detailed microscopic and 
macroscopic properties of different LZ systems. 

In the fundamental problems of chemical dynamics and molecular spectroscopy, the transitions from the initial to 
final states can be treated as a certain motion along the potential energy surfaces of the system under consideration. 
These surfaces in own turn are usually determined within the Born - Oppenheimer approximation (see section II). 
However, the approximation becomes inadequate for the excited vibrational states, when their energies are of the order 
of electronic inter level energy spacing or near the dissociation limit. In the both cases the non-adiabatic transitions 
should be taken into account, and the most of the non-radiative processes occur owing to this non-adiabaticity. The 
typical examples investigated in the monography [80], are so-called pre-dissociation, singlet-triplet or singlet-singlet 
conversion, and vibrational relaxation phenomena. 

Slow atomic collisions provide other examples of the non-adiabatic transitions between electronic states, where the 
time dependence of the states is determined by the distance and by the relative velocity of the colliding particles [33] . 
Some examples of the non-adiabatic transitions relevant for semiconductor physics can be found in [81], for nuclear or 
elementary particle physics in [82] and for laser or non- linear optic physics in [83] - [86]. The latter topic of course is 
of interest in its own right but also as an illustration of novel and fundamental quantum effects related to LZ model. 
The off-diagonal electronic state interactions are arisen for this case from the dipole forces. For relatively short laser 
pulses, it leads to the time dependent LZ problem for two electronic states, detailed considered in our paper (see also 
the laser optic formulation in [83] - [85] ) . The probability to find the system in the upper state after a single resonant 
passage, can be computed in the frame work of the LZ model. The latter point is related to one important aspect of 
the LZ problem, namely dissipative and noisy environments. When external actions (say fields) driving LZ transitions 
are reversed from large negative to large positive values, the dissipation reduces tunneling, that is the system remains 
in the ground state, or by other words, the thermal excitation from the ground state to the excited one, suppresses 
such adiabatic transitions. However, in the case of the field swept from the resonance point, the tunneling probability 
becomes larger in the presence of the dissipation (see e.g., [67]). The increasing precision of experimental tests in 
the femtosecond laser pulse range enables to excite well defined molecular states and to study their evolution in time 
using the second probing laser beam [17]. 
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Figure Captions. 

Fig. 1 

Stokes (dashed) and anti-Stokes (solid) lines for a pair of close linear turning points replaced by one second order 
turning point; (a) - classically forbidden region; (b) classically accessible region. 
Fig. 2 

Adiabatic (3 , 4) and diabatic (1 , 2) potentials for LZ problem. 
Fig. 3 

Stokes (dashed) and anti-Stokes (solid) lines in the vicinity of: (a) conjugated bifurcation points iirc; (b) diabatic 
potentials crossing point X = 0. 
Fig. 4 

Stokes (dashed) and anti-Stokes (solid) lines for linear turning points corresponding classically forbidden (a), inter- 
mediate (b), and accessible (c) energy regions of LZ problem. 
Fig. 5 

Relative placement of the adiabatic levels; (a) U12 > U^2 ; (b) U12 < U^2 . (^^12 = {^/2){n^F'^ /AmY/^). 
Fig. 6 

E, U12 phase diagram (/ - tunneling region, II - over-barrier region, and two intermediate energy regions /// and 
///' are separated by the line v* = 0.325. 
Fig. 7 

Transition matrix element Mn as a function of U-12/U12, computed at a = 0: on the boundary between tunneling 
and intermediate energy regions (a); a,t E = (b); on the boundary between the intermediate and over-barrier 
regions (c); 

lines 1,2,3, l',2',3' , 1",2",3" computed for corresponding energy regions using (5.19), (5.25), and (5.36), respec- 
tively. 
Fig. 8 

T versus E dependence for: (a) U12 = U^^ \ (b) U12 = 0.5C/i*2 ; (c) U12 = O.2WI2 ; 

stars mark the region III' boundaries, thin lines show results for the over-barrier and tunneling regions, bold lines 
for the intermediate energy region. 
Fig. 9 

Transmitted wave phase as a function of E in the over-barrier region at U\2 = U^2- 
Fig. 10 

Amplitudes of the decaying solutions at X > versus E for: (1) U12 = t^i*2 ; (2) U12 = 0.5Ui2 ; (3) U12 = 0.25Ui2- 
Fig. 11 

The diabatic level crossing phenomena: (a) crossing region; (b) bound initial and decay final states; (c) bound 
initial and final states. 
Fig. 12 

Level displacements versus U12 for two diabatic crossing potentials (l±X)^/2. Dashed lines show the intermediate 
energy region; dotted - dashed lines show displacements for the top and for the bottom of the adiabatic potentials. 
k, n, and n' are quantum numbers for the diabatic, and lower and upper adiabatic potentials. 

Fig. 13 

r„ versus U12 for the quasi stationary states at the diabatic potentials (1 + X)^/2 and (1/2) — X crossing; (a) 1 
- 4 are the level energies 0.042 , 0.125 , 0.208 , and 0.292 for the lower adiabatic potential; (b) 1' — 3' are the level 
energies 0.625 ; 0.708 ; 0792 for the upper adiabatic potential. 

Fig. 14 

Survival probability for the localized n = state; (a) b = 1500 , dashed lines U12 = 0.15 ; solid lines U12 = 0.21; 
(b) b = 1500, dashed lines U12 = 0.28 ; sohd lines U12 = 0.21. 
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